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1 Introduction 

In the previous lecture, we wrote the EFIE for an incident TE plane wave on a PEC surface. 
The solution was then obtained by some types of "intuitive" arguments, such as dividing the 
integration domain into small elements and supposing that the unknown does not vary too much 
over each elementary cell. 

We shall now see more rigorously what we actually did, and show that it was in fact a simple 
version of the Method of Moments. 

R. F. Harrington was the first to use the method of moments (MoM) in electromagnetics 
and his book remains a fundamental reference (and very easy to read!): 

R. F. Harrington, "Field Computation by Moment Method^ (is now available from IEEE Press). 



The situation we studied last time is depicted in Fig. 1. The integral equation (EFIE) we 
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Figure 1: Rough surface S separating two media. 



eventually obtained was: 
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Section 2. PEC surface with TE incident wave: EFIE 



where 

[n ■ V$(r-' )] U'=/(.o = U{x') , (3a) 

$i„c(x,/(x)) = 6(x), (3b) 
K{x,x) = g{xj{x),x\f{x)) , (3c) 

where U {x') is the unknown we are solving for. 

The important step we did from this integral, although it probably appeared straightforward, 
was to say that x can only take discrete values on the surface, thus defining TV intervals of length 
Ax: 

xe{xi}, i = l,...,N, (4) 

The assumptions were therefore: 

1. X G {xi} , i = 1, . . . , N. 

2. U{x') is constant on each interval. 

which eventually yielded the following system of equations (supposing that the problem related 
to the singularity of the Green's function has been accounted for): 

N 

^ ^ ^mn = hra ■ (5) 
n=l 

This is a matrix equation with the two indices m and n corresponding to: 
m: observation point unprimed coordinates. 

n: source point primed coordinates. 
Physical interpretation: 

• element (m, n) represents the effect of cell n on cell m. 

• element (m, m) represents the self-term. 

Mathematically: 

The steps we had to perform to from Eq. (2) to Eq. (5) are 

1. Write the unknown as 

U{x') = Y,Un5{Xn)^X, (6) 

n 

stating that now the unknowns become {Un} which are the amplitudes of the function. 
The integral equation becomes: 

r-L/2 

dx'K{x, x')U{x') = Ax V UnK{x, Xn) = b{x) . (7) 

n 
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2. Dot-multiply both sides by S{xm)' 



(8) 



n 



These two steps are at the basis of the method of moments: 

1. First step can actually be decomposed into two steps: 

(a) Mesh the structure (i.e. choose the intervals over which Un will be defined). 

(b) Expand the unknown U{x) into basis functions. 

2. The second steps concerns the observation: dot-multiply both sides of the equation by a 
TEST FUNCTION (or WEIGHTING FUNCTION). 

In the previous example: 

• basis functions: pulse basis function. 

• testing functions: we point-patch the integral equation at x = Xm (the method is therefore 
called point matching). 

This is a very simple, yet very widely used version of the method of moments. 

3 General considerations on MoM 

Let us consider the inhomogeneous equation: 



where £ is a linear operator, g is known, and / is to be determined. We shall now perform the 
two essential steps we have highlighted above. 

1. Let / be expanded in a series of functions: 



where q;„ are constant. The set /„ is called expansion function, or basis functions. 

Note that for an exact solution, the summation should be taken to oo, but has to be 
truncated in practice. 

2. It is assumed that a suitable inner product has been defined for the problem. Now, we 
define a set of weighting functions, or testing functions, wi , W2 , ■ ■ . , wn in the range of 
C, and take the inner product of the previous equation with Wm- 



A/) = 5, 



(9) 




(10) 



n 




(11) 



n 



4 



Section 4. A simple example for electrostatic 



The system can now be written in matrix form as: 

[Amn] [an] = [S'm] , 

where 



(12) 
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If the matrix [^rnn] is not singular, the unknowns are simply given by: 

[an] = [Amn]''^ [9m] , 



(14) 



and the original function / can be reconstructed using Eq. (10). We can now generalize the 
following definitions: 



The basis functions used previously are defined as: 
Pulse basis functions: /„ = < 



1 if X belongs to the interval n 
otherwise 



(15) 



• The testing (or weighting functions): 

Point matching = taking Dirac 6 functions as testing functions. 

4 A simple example for electrostatic 

The example is taken from the reference mentioned at the beginning of this document. 

Let us consider a square plate of side 2a lying on the z = plane with its center at the 
origin (see Fig. 2. Let a{x,y) represent the surface density on the plate, assumed to have zero 
thickness. The electrostatic potential $ at any point in space is given by 



2a 




Figure 2: Discretized square plate. 
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/a pa 
dx' I dy a{x' ,y')g{f,r) , 
-a J —a 

R=^/{x- x'Y + {y- y'Y + . 



(16a) 

(16b) 
(16c) 



The boundary condition \s ^ = V = constant on the plate. The integral equation for the 
problem is therefore 

a{x\y') 



/a ra 



(17) 



47re ^(^x-x'Y + {y-y'Y + z'^ ' 
where the unknown to determine is a(x',y'). Let us perform the three steps mentioned before: 
1. Mesh the structure: divide the plate into N squares of size 26 (see Fig. 2). 



2. Basis functions: let us choose 



N 



a{x', y') ~ ^ a„ /„(x', y') with /„(x', y') 



1 on ASn 

on AS"^, m^n. 



(18) 



3. Test functions: we choose to satisfy the integral equation at the mid-point (x^, y^) of 
each lS.Syn'- 

Wm = S{x-Xm)S{y-ym)- (19) 



With these three steps, we construct the matrix as (^; = 0): 



[Amn]= / dx' I 



dy' 



(20) 



lAxr, J Ayr, 47re ^(xm - x')2 + (y^ - y')^ 

It is obvious to see that this integral is singular at {xm,ym) S ASm- In this simple case 
fortunately, we can perform the integration analytically (this is not always the case), and write: 

dx' / dy' , = — ln(l + V2) 

AS„. 
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2^ _ 

AneRmn i^t^ix^ - x'f + (y^ - y'f ' 
To rewrite this with the language of linear space: 

j{x,y) = a{x,y) , 

g{x, y) = V on the plate (the discretization gives: g„ 

, 1 f{x',y') 



m ^ n . 



(21a) 
(21b) 



) 



(22a) 
(22b) 

, (22c) 

Aire y (x - x'y + {y - y'f 

Note: if we add another plate under the existing one, at z = —2c/, with a potential — we 
build a new problem that can be analyzed in two ways: 
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Section 5. Vectorial MoM 



1. By meshing both plates (ie. meshing everything, which is of course always possible). This 
will yield matrix twice as big as the previous one, to solve for twice as many unknowns. 

2. By using the image theory, and saying that the new problem is equivalent to the one of a 
unique plate on top of a ground plane at z = —d. In that case, we only have to change 
the Green's function to take the ground plane into account, and we keep the same number 
of unknowns as in the initial problem. When possible, this solution is better because 
computationally less expensive (analytically more expensive). 

Basically, a general trend is to have a Green's function that represents as much as possible 
the environment and to mesh only those parts that are external to the environment. This is in 
fact the reason why people are looking for Green's functions in layered media, periodic media, 
etc. 



5 Vectorial MoM 

We can of course apply the MoM to the vectorial case, like for example the equation: 

E{r) = icofi [ ds' ^e{f,f') ■ J{f') . (23) 
Js 

The general expansion of the current will be: 

J(f') = 5]a„/„(f'), (24) 

n 

yielding 

E{f) = iufi ^CXn [ ds' ^e(r, f') • fn{r') . (25) 

The third step is to dot-multiply the equation with a testing function hjn and integrate over 
the cell surface {i.e. perform the inner product): 



/ E{r) ■ hm{r) = iu/i'^ an ds' / ds hm{r) ■ Ge{r,r') ■ fn{r' 



) • (26) 



The double double-integral on the right-hand side of this equation is known as the 
"impedance term" since we can cast this system of equation into a matrix representation as: 

[Em] = [Zmn] [^n] • (27) 



6 Other basis and testing functions 

The advantage of the MoM over purely numerical methods is that there is still a large part that 
remains analytic (like the Green's functions for example). Yet, it remains a numerical method 
based on a matrix inversion technique and therefore, convergence issues need to be examined. 

The convergence of the MoM is closely related to the choice of basis functions and, although 
to a lesser extend, to the choice of testing functions. There are essentially two families of basis 
functions: 
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1. Entire domain basis functions: using these functions to expand the unknowns is analogous 
to a Fourier expansion or to a modal expansion. These types of functions yield a good 
convergence of the method but are not versatile since the geometry need be regular in order 
to have the modes defined. Note that in this case there is no use to mesh the geometry. 

2. Sub-domain basis functions: they rely on a proper meshing of the geometry, which can be 
rectangular, triangular, etc. The choice of basis functions is here very wide, from Dirac S 
(like for the weighting functions shown in this document), pulses (basis functions shown 
in this document), piecewise linear, etc. 

Finally, we can mention that point matching, which is easy to grasp and straightforward to 
implement, may not yield an optimal convergence. In most of the applications, the Galerkin 
technique is better, which consists in choosing the same testing functions as the basis functions. 
This applies to both sub-domain and entire domain functions. 



